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Abstract. The tight-binding model of electrons in graphene is reviewed. We derive 
low-energy Hamiltonians supporting massless Dirac-like chiral fermions and massive 
chiral fermions in monolayer and bilayer graphene, respectively, and we describe 
how their chirality is manifest in the sequencing of plateaus observed in the integer 
quantum Hall effect. The opening of a tuneable band gap in bilayer graphene in 
response to a transverse electric field is described, and we explain how Hartree 
theory may be used to develop a simple analytical model of screening. 



1.1 Introduction 

More than sixty years ago, Wallace [I] modeled the electronic band structure 
of graphene. Research into graphene was stimulated by interest in the proper- 
ties of bulk graphite because, from a theoretical point of view, two-dimensional 
graphene serves as a building block for the three-dimensional material. Fol- 
lowing further work, the tight-binding model of electrons in graphite, that 
takes into account coupling between layers, became known as the Slonczewski- 
Weiss-McClure model [21 EH]. As well as serving as the basis for models 
of carbon-based materials including graphite, buckyballs, and carbon nan- 
otubes [SliniElElElIinilH] , the honeycomb lattice of graphene has been used 
theoretically to study Dirac fermions in a condensed matter system p~2ll 13] . 
Since the experimental isolation of individual graphene flakes [13] . and the 
observation of the integer quantum Hall effect in monolayers [15(116] and bi- 
layers [17) . there has been an explosion of interest in the behavior of chiral 
electrons in graphene. 

This Chapter begins in Sect. ll.2l with a description of the crystal structure 
of monolayer graphene. Section [1.31 briefly reviews the tight-binding model of 
electrons in condensed matter materials [THllllj . and Sect. ll.4l describes its ap- 
plication to monolayer graphene [ITJ[l9l[2Q] . Then, in Section II. 5\ we explain 
how a Dirac-like Hamiltonian describing massless chiral fermions emerges from 
the tight-binding model at low energy. The tight-binding model is applied to 
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bilayer graphene in Sect. 11.61 and Sect. 11.71 describes how low-energy elec- 
trons in bilayers behave as massive chiral quasiparticles [T7][2T]. In Sect. 11.81 
we describe how the chiral Hamiltonians of monolayer and bilayer graphene 
corresponding to Berry's phase ir and 2n, respectively, have associated four- 
and eight-fold degenerate zero-energy Landau levels, leading to an unusual 
sequence of plateaus in the integer quantum Hall effect [15, 161117] . 

Section 11.91 discusses an additional contribution to the low-energy Hamil- 
tonians of monolayer and bilayer graphene, known as trigonal warping [4]I22| 
H21IM1I2S1H1I2I] > that produces a Liftshitz transition in the band structure of 
bilayer graphene at low energy. Finally, Sect. ITTTOl describes how an external 
transverse electric field applied to bilayer graphene, due to doping or gates, 
may open a band gap that can be tuned between zero up to the value of the 
interlayer coupling, around three to four hundred meV [3TJ|2BJ|27]- Hartree 
theory and the tight-binding model are used to develop a simple model of 
screening by electrons in bilayer graphene in order to calculate the density 
dependence of the band gap [28] . 



1.2 The crystal structure of monolayer graphene 
1.2.1 The real space structure 

Monolayer graphene consists of carbon atoms arranged with a two-dimensional 
honeycomb crystal structure as shown in Fig. ll.lf a). The honeycomb structure 
[TSJ[TT] consists of the hexagonal Bravais lattice, Fig. Il.lf b). with a basis of 
two atoms, labeled A and B, at each lattice point. 

Throughout this Chapter, we use a Cartesian coordinate system with x 
and y axes in the plane of the graphene crystal, and a z axis perpendicular 
to the graphene plane. Two-dimensional vectors in the same plane as the 
graphene are expressed solely in terms of their x and y coordinates, so that, 
for example, the primitive lattice vectors of the hexagonal Bravais lattice, 
Fig. Il.lf b) , are ai and a 2 where 

and a = |ai| = |&2 1 is the lattice constant. In graphene, a = 2.46 A [TTj . 
The lattice constant is the distance between unit cells, whereas the distance 
between carbon atoms is the carbon-carbon bond length acc — a/V3 = 
1.42 A. Note that the honeycomb structure is not a Bravais lattice because 
atomic positions A and B are not equivalent: it is not possible to connect 
them with a lattice vector R = mai + 712^2 where n\ and ri2 are integers. 
Taken alone, the A atomic positions (or, the B atomic positions) make up an 
hexagonal Bravais lattice and, in the following, we will often refer to them as 
the l A sublattice' (or, the 'B sublattice'). 
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Fig. 1.1. (a) The honeycomb crystal structure of monolayer graphene where white 
(black) circles indicate carbon atoms on A (B) sites and straight lines indicate a 
bonds between them. Vectors ai and a2 are primitive lattice vectors of length equal 
to the lattice constant a. The shaded rhombus is a unit cell containing two atoms, 
one A and one B. (b) Crosses indicate lattice points of the hexagonal Bravais lattice. 
The honeycomb structure in (a) consists of the hexagonal Bravais lattice [shown in 
(b)] with a basis of two atoms, one A and one B, at each lattice point. 



1.2.2 The reciprocal lattice of graphene 

Primitive reciprocal lattice vectors t>i and D2 satisfying aibi = a2b2 = 2ir 
and aib2 = a2bi = are given by 

/2?r 2tt \ /2vr 2tt \ . . 

b 1= — ,— , b 2 = — ,--=- . (1.2) 



a ' y/3aj ' V a ' ^ a 

The resulting reciprocal lattice is shown in Fig. ll.2[ which is an hexagonal Bra- 
vais lattice. The first Brillouin zone is hexagonal, as indicated by the shaded 
region in Fig. 11.21 

1.2.3 The atomic orbitals of graphene 

Each carbon atom has six electrons, of which two are core electrons and four 
are valence electrons. The latter occupy 2s, 2p Xl 2p y , and 2p z orbitals. In 
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Fig. 1.2. The reciprocal lattice of monolayer graphene where crosses indicate recip- 
rocal lattice points, and vectors bi and b2 are primitive lattice vectors. The shaded 
hexagon indicates the first Brillouin zone. 

graphene, the orbitals are sp 2 hybridized, meaning that two of the 2p orbitals, 
the 2p x and 2p y that lie in the graphene plane, mix with the 2s orbital to 
form three sp 2 hybrid orbitals per atom, each lying in the graphene plane 
and oriented 120° to each other [11]. They form a bonds with other atoms, 
shown as straight lines in the honeycomb crystal structure, Fig. 11.1( a). The 
remaining 2p z orbital for each atom lies perpendicular to the plane, and, 
when combined with the 2p z orbitals on adjacent atoms in graphene, forms a 
7T orbital. Electronic states close to the Fermi level in graphene are described 
well by a model taking into account only the it orbital, meaning that the 
tight-binding model can include only one electron per atomic site, in a 2p z 
orbital. 



1.3 The tight-binding model 

We begin by presenting a general description of the tight-binding model for 
a system with n atomic orbitals <f>j in the unit cell, labeled by index j = 
1 . . . n. Further details may be found in the book by Saito, Dresselhaus, and 
Dresselhaus |llj . It is assumed that the system has translational invariance. 
Then, the model may be written using n different Bloch functions <Pj (k, r) 
that depend on the position vector r and wave vector k. They are given by 



1 N 

# i (k J r) = 7 =XI eik " R ' , ^( r - R A0 - 

V i=l 



(1.3) 
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where the sum is over N different unit cells, labeled by index i = 1 . . . N, and 
Rj i denotes the position of the jth orbital in the ith unit cell. 

In general, an electronic wave function (k, r) is given by a linear super- 
position of the n different Bloch functions, 



^(k,r)=^c jV (k)^(k,r), (1.4) 
i=i 



where Cj t i are coefficients of the expansion. The energy Ej(k) of the jth band 
is given by 



EM = (L5) 



where H is the Hamiltonian. Substituting the expansion of the wave function 
(|1.4|) into the energy gives 



E"; HtiCjjCjl 



(1.7) 



where transfer integral matrix elements Hn and overlap integral matrix ele- 
ments Su are defined by 

H a = ($i\H\$i) , S a = ($i\$i). (1.8) 

We minimize the energy Ej with respect to the coefficient c* m by calculating 
the derivative, 

dEj _ Er ' H m iCji _ E"; H ii c *^ c ji E" gmfCjj 



The second term contains a factor equal to the energy Ej itself, (|1.7p . Then, 
setting dEj/dc* m — and omitting the common factor Suc^Cji gives 

n n 
^ HmlCjl = Ej ^ SmlCjl ■ (1-10) 
{=1 {=1 

This can be written as a matrix equation. Consider the specific example of 
two orbitals per unit cell, n = 2. Then, we can select the possible values of m 
(cither m = 1 or m = 2) and write out the summation in (| 1 . 10[) explicitly: 

m = 1 HnCj! + H 12 c j2 = Ej (SuCji + S 12 c o2 ) , (1-H) 

m = 2 H 2 iCji + H 2 2C j2 = Ej (S 2 ic jl + S 22 c j2 ) . (1.12) 
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These two equations may be combined into a matrix equation, 

5*11 5*12 



#11 H12 

H21 H22 



Cjl 
Cj 2 



Ei 



S21 S22 



Cjl 
C 3 2 



(1.13) 



For general values of n, defining H as the transfer integral matrix, S as the 
overlap integral matrix and f/'j as a column vector, 



H = 



I Hll H12 ■ ■ ■ Hi n \ 
H2I H22 ■ ■ ■ H2n 



S 



\ Snl S n 



5*12 • 


■ Sin ^ 




( C ^ \ 


■5*22 ' 


■ S211 




C j2 






> $3 = 




S n 2 ■ 


Snn J 




\ C 3 n / 



\ H n l H n 2 ■ ■ ■ H n n ) 

allows the relation (|1 . 10|) to be expressed as 

Hipj = EjSipj ■ 

The energies Ej may be determined by solving the secular equation 

det {H-EjS) =0, 



(1.14) 



(1.15) 



(1.16) 



once the transfer integral matrix H and the overlap integral matrix S are 
known. Here, 'det' stands for the determinant of the matrix. In the following, 
we will omit the subscript j — 1 ... n in (|1.15j) , (I1.16|) , bearing in mind that 
the number of solutions is equal to the number of different atomic orbitals per 
unit cell. 



1.4 The tight-binding model of monolayer graphene 

We apply the tight-binding model described in Sect. ll.3l to monolayer graphene, 
taking into account one 2p z orbital per atomic site. As there are two atoms in 
the unit cell of graphene, labeled A and B in Fig. 11.11 the model includes two 
Bloch functions, n — 2. For simplicity, we replace index j = 1 with j = A, and 
j = 2 with j = B. Now we proceed to determine the transfer integral matrix 
H and the overlap integral matrix S. 

1.4.1 Diagonal matrix elements 

Substituting the expression for the Bloch function (11.31) into the definition 
of the transfer integral (jl.8|) allows us to write the diagonal matrix element 
corresponding to the A sublattice as 

N N 

H ^ = Jj EE e^ R ^- RA ^(4>A (r - R^) \H\<j>A (r - R^)), (1.17) 

i=i j=i 



1 Electronic properties of monolayer and bilayer graphene 



7 



where k = (k x ,k v ) is the wave vector in the graphene plane. Equation (|1.17[) 
includes a double summation over all the A sites of the lattice. If we assume 
that the dominant contribution arises from the same site j = i within every 
unit cell, then: 

1 N 

Ha a « jz Y,^a (r - Ri,j) (r - R A ,i)) , (1.18) 

i=i 

The matrix element {4> a\H\4> a) within the summation has the same value on 
every A site, i.e. it is independent of the site index i. We set it to be equal to 
a parameter 

e 2p = {<t>A (r - R A ,i) \H\cj>A (r - R A ,i)) , (1-19) 

that is equal to the energy of the 2p 2 orbital. Then, keeping only the same 
site contribution, 



1 N 

i=l 

It is possible to take into account the contribution of other terms in the double 
summation (jl.171) . such as next-nearest neighbor contributions [29ll30j . They 
generally have a small effect on the electronic band structure and will not be 
discussed here. The B sublattice has the same structure as the A sublattice, 
and the carbon atoms on the two sublattices are chemically identical. This 
means that the diagonal transfer integral matrix element corresponding to the 
B sublattice has the same value as that of the A sublattice: 

H BB = H AA ^e 2p . (1.21) 

A calculation of the diagonal elements of the overlap integral matrix pro- 
ceeds in a similar way as for those of the transfer integral. In this case, the 
overlap between a 2p z orbital on the same atom is equal to unity, 

(<t> A (r - R Aii ) \<j> A (r - R Aii )) = 1 . (1.22) 

Then, assuming that the same site contribution dominates, 

N N 

X 



S ^ = JfT,T, e lk ' (R ^- R ^ } (4>a (r - R A ,i) \4>a (r - Raj)) , (1.23) 

i=i j=i 

1 N 

- ]T(<^4 (r - R A ,i) \4>a (r - R A ,i)) , (1-24) 

i=l 

^ N 

mT, 1 (1-25) 



N ■ 

i=l 
N 

.V ■ 
1=1 

1. (1.26) 
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Fig. 1.3. The honeycomb crystal structure of monolayer graphene. In the nearest- 
neighbor approximation, we consider hopping from an A site (white) to three adja- 
cent B sites (black), labeled Bi, B2, B3, with position vectors Si, S2, S3, respectively, 
relative to the A site. 

Again, as the B sublattice has the same structure as the A sublattice, 

Sbb =S aa = 1. (1.27) 



1.4.2 Off-diagonal matrix elements 

Substituting the expression for the Bloch function (jl.3p into the definition of 
the transfer integral (jl.8p allows us to write an off-diagonal matrix element as 

1 N N 

H ^ b= nJ2Y1 e*- (RB ^ RAi) <^ (r - R A ,i) (r - R Bj3 -)>- (1-28) 

t=l j=X 

It describes processes of hopping between the A and B sublattices, and con- 
tains a summation over all the A sites (i — 1 . . . N) at positions Ra,« and all 
the B sites (j = 1 . . . N) at R_bj. 

In the following, we assume that the dominant contribution to the off- 
diagonal matrix element (|1.28j) arises from hopping between nearest neighbors 
only. If we focus on an individual A atom, i.e. we consider a fixed value of the 
index i, we see that it has three neighboring B atoms, Fig. 11.31 that we will 
label with a new index I (I = 1 ... 3). Each A atom has three such neighbors, 
so it is possible to write the nearest-neighbors contribution to the off-diagonal 
matrix element (|1 .28|) as 

1 N 3 

Hab ~ n £ 2 e i;k - (RB ' ! - R ^- ) (0A (r - Ra,0 (r - R B| ,)) ■ (1-29) 

i=l 1=1 
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The matrix element between neighboring atoms, {4>aO~L\4>b) , has the same 
value for each neighboring pair, i.e. it is independent of indices i and I. We 
set it equal to a parameter, t — (<j>A (r — R. A ,i) \H\4>b ( r — Rb,0)- Since t 
is negative it is common practice to express it in terms of a positive 
parameter 70 = — t, where 

70 = -{fa (r - K A ,i) \U\ct> B (r - R B ,i)) ■ (1-30) 
Then, we write the off-diagonal transfer integral matrix element as 

1 n 3 

HAB^—^^^'-^ho, (1.31) 
j=l 1=1 

N 3 

= -|EE elk '^-^(k), (1-32) 

i=l i=l 

/(k) = ^e ik -^, (1.33) 
i=i 

where the position vector of atom Bi relative to the Ai atom is denoted 
Si = ~R.b,i ~ Rvi.i, and we used the fact that the summation over the three 
neighboring B atoms is the same for all Ai atoms. 

For the three B atoms shown in Fig. 11.31 the three vectors are 

( i \ » fa a \ „ fa a \ 

Note that |<$i| = |<$2| = \8a\ = a/V3 is the carbon-carbon bond length. Then, 
the function / (k) describing nearest-neighbor hopping may be evaluated as 

/(k) = ]ry k ^, (1.35) 

1=1 

_ e ik v a/y/3 _|_ e ife x a/2 e -ifc H a/2v / 3 ^_ e ~ik x a/2 e -ifc H a/2\/3 M gg-j 

= e ik y a/V3 + 2e^ fc « a / 2 ^cos (k x a/2) . (1.37) 
The other off-diagonal matrix element Hba is the complex conjugate of Hab' 
Hab « -70/ (k) , H BA « -70/* (k) . (1.38) 

A calculation of an off-diagonal element of the overlap integral matrix 
proceeds in a similar way as for the transfer integral: 



S *b = i^^e ik -( R ^- R ^)(^(r-R^)|^(r-R BJ )), (1-39) 

i=i j=i 

1 N 3 

~ «EE e*-^-'- 11 ^^^ (r - Rx,i) \4>b (r - R B ,i)> , (1-40) 



N ■ 

i=i i=i 



so/(k), (1.41) 
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where the parameter so = (4>a ( r — R-A,i) \4>b ( r — R-B,i))' an< ^ = ^ab = 
s f* (k). The presence of non-zero Sq takes into account the possibility that 
orbitals on adjacent atomic sites are not strictly orthogonal. 



1.4.3 The low-energy electronic bands of monolayer graphene 

Summarizing the results of this section, the transfer integral matrix elements 
(|1.2ip and (|1.38|) . and the overlap integral matrix elements (|1 .27|) and (|1.41|) 
give 

Hx=( £ S^" 70/(k) V S 1 =( L^ S0f , {k) ) , (1-42) 
V~7o/*(k) e 2p J \s f*(k) 1 J 

where we use the subscript ' 1 ' to stress that these matrices apply to monolayer 
graphene. The corresponding energy E may be determined by solving the 
secular equation det (Hi — ESi) — 0, (|1.16[) : 

^- ( 7o + Es ) f* (k) e 2p -E J v ' 

(E-e 2p ) 2 -([E-e 2p ]s a + e 2pSo + lQ ) 2 \f{k)\ 2 =0. (1.44) 

Solving this quadratic equation yields the energy: 

_ e 2p ± 70 |/(k)| 
E± ~ lT*o|/(k)| ■ (L45) 

This expression appears in Saito et al [llj . where parameter values 70 = 
3.033 eV, so = 0.129, e 2p — are quoted. The latter value (e 2p — 0) means 
that the zero of energy is set to be equal to the energy of the 2p z orbital. The 
resulting band structure E± is shown in Fig. ll.4l in the vicinity of the Brillouin 
zone. A particular cut through the band structure is shown in Fig. 11.51 where 
the bands are plotted as a function of wave vector component k x along the 
line k y = 0, a line that passes through the center of the Brillouin zone, labeled 
r, and two corners of the Brillouin zone, labeled K + and K_ (see the inset 
of Fig. II .5[) . The Fermi level in pristine graphene is located at zero energy. 
There are two energy bands, that we refer to as the conduction band (£+) 
and the valence band (-E-). The interesting feature of the band structure is 
that there is no band gap between the conduction and valence bands. Instead 
the bands cross at the six corners of the Brillouin zone, Fig. 11.41 The corners 
of the Brillouin zone are known as K points, and two of them are explicitly 
labeled K + and K_ in Fig. 11.41 Near these points, the dispersion is linear 
and electronic properties may be described by a Dirac-like Hamiltonian. This 
will be explored in more detail in the next section. Note also that the band 
structure displays a large asymmetry between the conduction and valence 
bands that is most pronounced in the vicinity of the F point. This arises from 
the non-zero overlap parameter s appearing in (|1.45[) . 




Fig. 1.4. The low-energy band structure of monolayer graphene Eq. (|1.45p tak- 
ing into account nearest-neighbor hopping with parameter 70 = 3.033 eV, nearest- 
neighbor overlap parameter so = 0.129, and orbital energy e2 P = [IT]. The plot 
shows the bands calculated in the vicinity of the first Brillouin zone, with conduc- 
tion and valence bands touching at six corners of the Brillouin zone, two of them 
are labeled K+ and K- . Label r indicates the center of the Brillouin zone. 

The tight-binding model described here cannot be used to determine the 
values of parameters such as 70 and sq. They must be determined either by 
an alternative theoretical method, such as density-functional theory, or by 
comparison of the tight-binding model with experiments. Note, however, that 
the main qualitative features described in this chapter do not depend on the 
precise values of the parameters quoted. 

1.5 Massless chiral quasiparticles in monolayer graphene 
1.5.1 The Dirac-like Hamiltonian 

As described in the previous section, the electronic band structure of mono- 
layer graphene, Figs. I1.4[ 11.51 is gapless, with crossing of the bands at points 
K + and K- located at corners of the Brillouin zone. In this section, we show 
that electronic properties near these points may be described by a Dirac-like 
Hamiltonian. 
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Fig. 1.5. The low-energy band structure of monolayer graphene Eq. (|1.45[) tak- 
ing into account nearest-neighbor hopping with parameter 70 = 3.033 eV, nearest- 
neighbor overlap parameter sq = 0.129, and orbital energy t2 P = [11] . The plot 
shows a cut through the band structure Fig. 11.41 plotted along the k x axis inter- 
secting points K_, r, and K + in the Brillouin zone, shown as the dotted line in the 
inset. 

Although the first Brillouin zone has six corners, only two of them are 
non-equivalent. In this Chapter, we choose points K + and K_, Figs. 11.41 11.51 
as a non-equivalent pair. It is possible to connect two of the other corners 
to K + using a reciprocal lattice vector (hence, the other two are equivalent 
to 2£+), and it is possible to connect the remaining two corners to using 
a reciprocal lattice vector (hence, the remaining two are equivalent to K-) 1 
but it is not possible to connect K+ and K~ with a reciprocal lattice vector. 
To distinguish between K+ and K-, we will use an index £ = ±1. Using the 
values of the primitive reciprocal lattice vectors bi and b2, (11.21) . it can be 
seen that the wave vector corresponding to point K% is given by 



Note that the K points are often called 'valleys' using nomenclature from 
semiconductor physics. 

In the tight-binding model, coupling between the A and B sublattices is 
described by the off-diagonal matrix element Hab, (ll-38[) , that is proportional 
to parameter 70 and the function /(k), (| 1 . 35|) . Exactly at the point, k = 
K^, the latter is equal to 




(1.46) 
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/ (Kg) = e° + e l ^' z + e - l ^l z = . (1.47) 

This indicates that there is no coupling between the A and B sublattices 
exactly at the Kg point. Since the two sublattices are both hexagonal Bravais 
lattices of carbon atoms, they support the same quantum states, leading to a 
degeneracy point in the spectrum at Figs. 11.41 11.51 

The exact cancelation of the three factors describing coupling between 
the A and B sublattices, (jl.47l) . no longer holds when the wave vector is not 
exactly equal to that of the point. We introduce a momentum p that is 
measured from the center of the point, 

p = ftk-ftK e . (1.48) 

Then, the coupling between the A and B sublattices is proportional to 

/ (k) = e *pWVS» + 2e -i Pya /2Vzn cos f??£ + m\ f (L49) 

~ — ~ l Py) > ( 1 - 51 ) 

where we kept only linear terms in the momentum p = (p x ,Py), an ap- 
proximation that is valid close to the point, i.e. for pa/fi <C 1, where 
p = |p| = (p^ +Py) 1 ^ 2 - Using this approximate expression for the function 
/ (k), the transfer integral matrix (|1.42p in the vicinity of point K$ becomes 

^"UW'o"") ■ (152) 

Here, we used ti v = [llj which defines the zero of the energy axis to coincide 
with the energy of the 2p z orbital. The parameters a and 70 were combined 
into a velocity v defined as v — VSajo / (2ft) • 

Within the linear- in-momentum approximation for / (k), (|1.51[) . the over- 
lap matrix Si may be regarded as a unit matrix, because its off-diagonal el- 
ements, proportional to so, only contribute quadratic- in- momentum terms to 
the energy E±, (jl.45j) . Since S± is approximately equal to a unit matrix, (|1.15p 
becomes Hiip = Eip, indicating that H±, (|1.52[) . is an effective Hamiltonian 
for monolayer graphene at low-energy. The energy eigenvalues and eigenstates 
of Hi are given by 



E± = ±vp , V± = — [ ^ v ) e lp - r/h , (1-53) 



y/2 {±& 

where ± refer to the conduction and valence bands, respectively. Here ip is 
the polar angle of the momentum in the graphene plane, p = (p x ,Py) = 
p (cos ip, sin tp) . 
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Fig. 1.6. Schematic representation of the pseudospin degree of freedom: (a) elec- 
tronic density solely on the A sublattice can be viewed as a pseudospin 'up' state, 
whereas (b) density solely on the B sublattice corresponds to a pseudospin 'down' 
state; (c) in graphene, electronic density is usually shared equally between A and B 
sublattices, so that the pseudospin part of the wave function is a linear combination 
of 'up' and 'down,' with amplitudes dependent on the direction of the electronic mo- 
mentum p; (d) at valley K+, the pseudospin (<r) e in the conduction band is parallel 
to the momentum, whereas the pseudospin (ct)h in the valence band is anti-parallel 
to the momentum. 

1.5.2 Pseudospin and chirality in graphene 

The effective Hamiltonian (|1.52l) and eigenstates (11.53)) in the vicinity of the 
point have two components, reminiscent of the components of spin-1/2. 
Referring back to the original definitions of the components of the column vec- 
tor ip, (|1.4[) and (I1.14j) . shows that this is not the physical spin of the electron, 
but a degree of freedom related to the relative amplitude of the Bloch function 
on the A or B sublattice. This degree of freedom is called pseudospin. If all 
the electronic density was located on the A sublattice, Fig. |1.6f a). this could 
be viewed as a pseudospin 'up' state (pointing upwards out of the graphene 
sheet) |t) = (1,0) T , whereas density solely on the B sublattice corresponds 
to a pseudospin 'down' state (pointing downwards out of the graphene sheet) 
||) = (0, 1) T , Fig. |1.6f b). In graphene, electronic density is usually shared 
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Fig. 1.7. Anisotropic scattering of chiral electrons in graphene: (a) angular depen- 
dence w((p) — cos 2 (ip/2) of the scattering probability off an A-B symmetric potential 
in monolayer graphene [9lllQII33j and (b) w(<p) = cos 2 (y>) in bilayer graphene [211137] . 



equally between A and B sublattices, Fig. 11.6( c). so that the pseudospin part 
of the wave function is a linear combination of 'up' and 'down,' and it lies in 
the plane of the graphene sheet. 

Not only do the electrons possess the pseudospin degree of freedom, but 
they are chiral, meaning that the orientation of the pseudospin is related to 
the direction of the electronic momentum p. This is reflected in the fact that 
the amplitudes on the A or B sublattice of the eigenstate (|1.53[) depend on 
the polar angle ip. It is convenient to use Pauli spin matrices in the A/B 
sublattice space, <Ji where i = 1 ... 3, to write the effective Hamiltonian (|1.52p 
as 



Hl,t = v {&xPx + VyPy) • (I- 54 ) 

If we define a pseudospin vector as er = (<j x , a y , a z ), and a unit vector as 
ni = (£ cos ip, sin ip, 0), then the Hamiltonian becomes H\£ = uper.rii, stress- 
ing that the pseudospin er is linked to the direction fix . The chiral operator 
er.ni projects the pseudospin onto the direction of quantization rii: eigen- 
states of the Hamiltonian are also eigenstates of cr.hi with eigenvalues ±1, 
tr.&iip± — ±ip±- An alternative way of expressing this chiral property of 
electrons is to explicitly calculate the expectation value of the pseudospin op- 
erator (er) = ((a x ), (a v ), (c z )) with respect to the eigenstate ip±, (|1.53l) . The 
result, (er) e y^ = ± (£ cost/?, sin 0), shows the link between pseudospin and 
momentum. For valley K + , the pseudospin in the conduction band (er) e is 
parallel to the momentum, whereas the pseudospin in the valence band (er)/, 
is anti-parallel to it, Fig. |1.6f d). 

If the electronic momentum p rotates by angle ip, then adiabatic evolution 
of the chiral wave function ip±, (11.531) . produces a matching rotation of the 
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vector n! by angle tp. For traversal of a closed contour in momentum space, 
corresponding to ip = 2tt, then the chiral wave function undergoes a phase 
change of it known as Berry's phase [3TJI32]- It can be thought of as arising 
from the rotation of the pseudospin degree of freedom. 

The chiral nature of low-energy electrons in graphene places an additional 
constraint on their scattering properties. If a given potential doesn't break the 
A-B symmetry, then it is unable to influence the pseudospin degree of free- 
dom which must, therefore, be conserved upon scattering. Considering only 
the pseudospin part of the chiral wave function ip± , (|1.53p , the probability to 
scatter in a direction tp, where tp = is the forwards direction, is proportional 
to w(tp) — \(ip±(tp)\tp±(0)}\ 2 ■ F° r monolayer graphene, w{tp) = cos 2 {tp/2), 
Fig. |1.7f a). This is anisotropic, and displays an absence of backscattering 
w(ir) = [SJ[Tni|33]: scattering into a state with opposite momentum is pro- 
hibited because it requires a reversal of the pseudospin. Such conservation of 
pseudospin is at the heart of anisotropic scattering at potential barriers in 
graphene monolayers 34,35 , known as Klein tunneling. 

1.6 The tight-binding model of bilayer graphene 

In this section, we describe the tight-binding model of bilayer graphene. To do 
so, we use the tight-binding model described in Sect. 1 1.31 in order to generalize 
the model for monolayer graphene discussed in Sect. 11.41 

We consider Bernal-stacked bilayer graphene [36,[T7l[21] (also called AB- 
stacked bilayer graphene). It consists of two parallel layers of carbon atoms, 
each arranged with a honeycomb arrangement as in a monolayer, that are 
coupled together, Fig. 11.61 There are four atoms in the unit cell, a pair Al, 
Bl, from the lower layer and a pair A2, B2, from the upper layer. In Bernal 
stacking, the layers are arranged so that two atoms, Bl and A2, are directly 
below or above each other, whereas the other two atoms, Al and B2, do not 
have a counterpart in the other layer. The primitive lattice vectors ai and a.2, 
and the lattice constant a are the same as for monolayer graphene, and the 
unit cell, shown in Fig. I1.6f a). has the same area in the x-y plane as in the 
monolayer. Therefore, the reciprocal lattice and first Brillouin zone are the 
same as in monolayer graphene, Fig. 11.21 The unit cell of bilayer graphene 
contains four atoms, and, if the tight-binding model includes one p z orbital 
per atomic site, there will be four bands near zero energy, instead of the two 
bands in monolayer graphene. 

Essential features of the low-energy electronic band structure may be de- 
scribed by a minimal tight-binding model including nearest-neighbor coupling 
7o between Al and Bl, and A2 and B2, atoms on each layer, and nearest- 
neighbor interlayer coupling 71 between Bl and A2 atoms that are directly 
below or above each other, 



71 = {<PA2 (r - TL A 2) \U\<t>BX (r - R B i)) 



(1.55) 
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Fig. 1.8. Schematic representation of the crystal structure of A_B-stacked bilayer 
graphene: (a) plan view with Al (white) and Bl atoms (black) on the lower layer, 
A2 (black) and B2 atoms (grey) on the upper layer. Vectors ai and &2 are primitive 
lattice vectors of length equal to the lattice constant a, and the shaded rhombus is a 
unit cell; (b) side view where the parameter 70 represents nearest-neighbor coupling 
within each layer, 71 nearest-neighbor coupling between the Bl and A2 atoms on 
different layers. 



Then, we can generalize the treatment of monolayer graphene, (ll.42|) . to write 
the transfer and overlap integral matrices of bilayer graphene, in a basis with 
components Al, Bl, A2, B2, as 



H = 



S = 



^2p 

-70/* (k) 



*o/* (k) 




7o/ (k) 

7i 


*o/(k) 
1 




V 





71 

-7o/* (k) 




so/(k) 






-7o/ (k) 

(-2], 



\ 







(1.56) 



(1.57) 



sof(k) 1 / 



The upper-left and lower-right 2x2 blocks describe behavior within the lower 
{AljBi) and upper (A2/B2) layers, respectively. The off-diagonal 2 x 2 blocks, 
containing parameter 71, describe interlayer coupling. 
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The band structure of bilayer graphene may be determined by solving the 
secular equation det (H — EjS) = 0, (I1.16[) . It is plotted in Fig. 11.91 for pa- 
rameter values 70 = 3.033eV, s — 0.129, e 2p = [TT] and interlayer coupling 
71 = 0.39eV. There are four energy bands, two conduction bands and two 
valence bands. Overall, the band structure is similar to that of monolayer 
graphene, Fig. 11.51 with each monolayer band split into two by an energy 
approximately equal to the interlayer coupling 71 |36j . The most interesting 
part of the band structure is in the vicinity of the K points [21] , as shown in 
the left inset of Fig. 11.91 which focuses in on the bands around K_. At the K 
point, one of the conduction (valence) bands is split away from zero energy by 
an amount equal to the interlayer coupling 71 (-71). The split bands originate 
from atomic sites Bl and A2 that have a counterpart atom directly above or 
below them on the other layer. Orbitals on these pairs of atoms (Bl and A2) 
are strongly coupled by the interlayer coupling 71 and they form a bonding 
and anti-bonding pair of bands, split away from zero energy. In the follow- 
ing, we refer to them as 'dimer' states, and atomic sites Bl and A2 are called 
'dimer' sites. The remaining two bands, one conduction and one valence band, 
touch at zero energy: as in the monolayer, there is no band gap between the 
conduction and valence bands. In the vicinity of the K points, the dispersion 
of the latter bands is quadratic E± oc ±|k — K^| 2 , and electronic properties of 
the low-energy bands may be described by an effective Hamiltonian describ- 
ing massive chiral particles. This will be explored in more detail in the next 
section. 



1.7 Massive chiral quasiparticles in bilayer graphene 
1.7.1 The low-energy bands of bilayer graphene 

To begin the description of the low-energy bands in bilayer graphene, we set 
so = 0, thus neglecting the non-orthogonality of orbitals that tends to become 
important at high energy. Then, the overlap matrix S, (|1.57j) . becomes a unit 
matrix, and H, (jl.56|) . is an effective Hamiltonian for the four bands of bilayer 
graphene at low-energy [2Tj : 



H 



( -70/ (k) \ 

-70/* (k) 7l 

71 -70/ (k) 

V -70/* (k) J 



(1.58) 



where we used e 2p = [TT to define the zero of the energy axis to coincide 
with the energy of the 2p z orbital. Eigenvalues of the Hamiltonian are given 
by 



E 



(a) 




±1. 



(1.59) 
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Fig. 1.9. The low-energy band structure of bilayer graphene taking into account 
nearest-neighbor hopping with parameter 70 = 3.033eV, nearest-neighbor overlap 
parameter so = 0.129, orbital energy £2 P = [11], and interlayer coupling 71 = 
0.39eV. The plot shows the bands calculated along the k x axis intersecting points 
K-, r, and K+ in the Brillouin zone, shown as the dotted line in the right inset. 
The left inset shows the band structure in the vicinity of the point K- . 

Over most of the Brillouin zone, where 47Q|/(k)| 2 ^> j 2 , the energy may be 
approximated as E± w ±(7o|/(k)|+cryi/2), meaning that the a = ±1 bands 
are approximately the same as the monolayer bands, (ll.45p . but they are split 
by the interlayer coupling 71. The eigenvalues E±\ (11.591) . describe two bands 
that are split away from zero energy by ±71 at the K point (where |/(k)| = 0) 
as shown in the left inset of Fig. 11.91 This is because the orbitals on the A2 
and B\ sites form a dimer that is coupled by interlayer hopping 71, resulting 
in a bonding and anti-bonding pair of states ±71. 

The remaining two bands are described by E± 1 '. Near to the point, 
pa/% <C 1, we replace the factor 7 |/(k)| with vp, (|1 .51[) : 



This formula interpolates between linear dispersion at large momenta (71 <C 
vp < 70) and quadratic dispersion E± ^ m ±v 2 p 2 /71 near zero energy where 
the bands touch. These bands arise from effective coupling between the or- 
bitals on sites, Al and B2, that don't have a counterpart in the other layer. 
In the absence of direct coupling between Al and B2, the effective coupling 
is achieved through a three stage process as indicated in Fig. 11.101 It can be 
viewed as an intralayer hopping between Al and Bl, followed by an interlayer 




(1.60) 
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Fig. 1.10. Schematic representation of the crystal structure of AB-stacked bilayer 
graphene illustrating the processes that contribute to effective coupling between Al 
(white) and B2 atoms (grey), in the presence of strongly-coupled 'dimer' sites Bl 
and A2 (black). The black arrowed line indicates the three stage process: intralayer 
hopping between Al and Bl, followed by an inter layer transition via the dimer sites 
Bl and A2, followed by another intralayer hopping between A2 and B2. 

transition via the dimer sites Bl and A2, followed by another intralayer hop- 
ping between A2 and B2. This effective coupling may be succinctly described 
by an effective low-energy Hamiltonian written in a two-component basis of 
p z orbitals on Al and B2 sites. 

1.7.2 The two-component Hamiltonian of bilayer graphene 

The effective two-component Hamiltonian may be derived from the four com- 
ponent Hamiltonian, (|1.58[) . using a Schrieffer- Wolff transformation [3"51l21| . 
In the present context, a straightforward way to do the transformation is to 
consider the eigenvalue equation for the four component Hamiltonian, (|1.58l) , 
as four simultaneous equations for the wave- function components cai, cbi, 

C-A2, CB2- 



Using the second and third equations, (II .62)) and (ll.63|) . it is possible to 
express the components on the dimer sites, cbi and CA2, in terms of the other 
two: 



Ecai + Jof (k) cbi 
7o/* (k) cai + Ec B i - 71CA2 

-7lCBl + Ec A 2 + 7o/ (k) C B 2 

70/* (k) C A 2 + Ec B 2 



0. 



0. 

0, 







(1.61) 
(1.62) 
(1.63) 
(1.64) 



cbi = 




CB2 + 



Ejpf* (k) 
lid 



cai, 



(1.65) 
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£70/ ft) , 70/* ft) f-i 

c -42 = 23 — C B2 H 3 — cai , (1.66) 

jfd 71 d 

where d = 1 — E 2 /^ 2 . Substituting these expressions into the first and fourth 
equations, (jl.6ip and (|1.64[) , produces two equations solely in terms of cai 
and CB2- Assuming \E\ -C I71I and (70/ (k) | -C |7i|, we use d ~ 1 and keep 
terms up to order I/71 only: 

Ec A1 + lof (k) CB2 = 0, (1.67) 
7i 

7o 2 (r ( k )) 2 CA1 + gC52 = . (1 . 68) 
71 



It is possible to express these two equations as a Schrodinger equation, H 2 ip 
Eip, with a t"w 
Hamiltonian 



Eip, with a two-component wave function ip = {cai,cb2) T an d two-component 



ff2 ' £ -~2^l(^ + ^) 2 )' (L69) 

where we used the approximation / (k) w — u (^p^ — ip^) /70, (11.51[) . valid for 
momentum pa/Ti <C 1 close to the JQ point, and parameters i; and 71 were 
combined into a mass m = 7i/(2w 2 ). 

The effective low-energy Hamiltonian of bilayer graphene, (jl.69p . resem- 
bles the Dirac-like Hamiltonian of monolayer graphene, (|1.52j) . but with a 
quadratic term on the off-diagonal instead of linear. The energy eigenvalues 
and eigenstates of H 2 are given by 

E ± =±^, ^ = -L( Te ^)e^A, (1.70) 

where ± refer to the conduction and valence bands, respectively. Here (p is 
the polar angle of the momentum in the graphene plane, p = (p x ,P y ) — 
p (cos <p, sin (f) . 



1.7.3 Pseudospin and chirality in bilayer graphene 

The two-component Hamiltonian (| 1 .69[) of bilayer graphene has a pseudospin 
degree of freedom [TTJET] related to the amplitude of the eigenstates (|1.70j) 
on the Al and B2 sublattice sites, where Al and B2 lie on different layers. 
If all the electronic density was located on the Al sublattice, Fig. II .lir a), 
this could be viewed as a pseudospin 'up' state (pointing upwards out of the 
graphene sheet) |f) = (1,0) T , whereas density solely on the B2 sublattice 
corresponds to a pseudospin 'down' state (pointing downwards out of the 
graphene sheet) ||) = (0, 1) T , Fig. ll.llT b). In bilayer graphene, electronic 
density is usually shared equally between the two sublattices, Fig. ll.llT c). so 
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Fig. 1.11. Schematic representation of the pseudospin degree of freedom in bilayer 
graphene: (a) electronic density solely on the Al sublattice on the lower layer can be 
viewed as a pseudospin 'up' state, whereas (b) density solely on the B2 sublattice on 
the upper layer corresponds to a pseudospin 'down' state; (c) in bilayer graphene, 
electronic density is usually shared equally between Al and B2 sublattices, so that 
the pseudospin part of the wave function is a linear combination of 'up' and 'down,' 
with amplitudes dependent on the direction of the electronic momentum p; (d) at 
valley K + , the pseudospin (cr) e in the conduction band is parallel to the quantization 
direction 112, whereas the pseudospin {(t)h in the valence band is anti-parallel to xi2- 
Direction £12 is related to the direction of momentum p, but turns in the x-y plane 
twice as quickly as it. 

that the pseudospin part of the wave function is a linear combination of 'up' 
and 'down,' and it lies in the plane of the graphene sheet. 

Electrons in bilayer graphene are chiral [TT1I2T| , meaning that the orienta- 
tion of the pseudospin is related to the direction of the electronic momentum 
p, but the chirality is different to that in monolayers. As before, we use Pauli 
spin matrices in the Al/ ' B2 sublattice space, Uj where i = 1 . . .3, to write the 
effective Hamiltonian (|1.69l) as 

#2,<; = -7^[°x (pI ~ Pl) + ZfrvPxPv] • (1-71) 

If we define a pseudospin vector as er = (a x , a y ,a z ), and a unit vector as 112 = 
— (cos2ip, £ sin 2(p,0), then the Hamiltonian becomes i?2,£ = (p 2 /2m) cr.n.2, 
stressing that the pseudospin er is linked to the direction n 2 . The chiral op- 
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erator <x.n 2 projects the pseudospin onto the direction of quantization n 2 : 
eigenstates of the Hamiltonian are also eigenstates of <x.n 2 with eigenvalues 
±1, (T.n2tp± — ±V J ±- I n bilayer graphene, the quantization axis n 2 is fixed 
to lie in the graphene plane, but it turns twice as quickly in the plane as the 
momentum p. If we calculate the expectation value of the pseudospin opera- 
tor (er) = ((<J X ) 7 (&y), {&z}) with respect to the eigenstate ip±, (|1.70j) . then the 
result (<r) e //j = =F (cos2<£>,£sin2</?,0), illustrates the link between pseudospin 
and momentum, Fig. ll.llT d). 

If the momentum p rotates by angle tp, adiabatic evolution of the chiral 
wave function tp±, (jl.70l) . produces a matching rotation of the quantization 
axis n 2 by angle 2tp, not tp as in the monolayer, Sect. 11.5.21 Thus, traversal 
around a closed contour in momentum space results in a Berry's phase |31II32| 
change of 2ir of the chiral wave function in bilayer graphene [TTIHI]- For 
Berry's phase 2-7T chiral electrons in bilayer graphene, (|1.70[) . the probability 
to scatter in a direction tp, where tp = is the forwards direction, is propor- 
tional to w(tp) = \(ip±((p)\ip±(0))\ 2 = cos 2 (ip) [2T1I37] as shown in Fig. OV) . 
This is anisotropic, but, unlike monolayers Fig. 11.7( a). does not display an 
absence of backscattering (w(tt) = 1 in bilayers): scattering into a state with 
opposite momentum is not prohibited because it doesn't require a reversal of 
the pseudospin. 

1.8 The integer quantum Hall effect in graphene 

When a perpendicular magnetic field is applied a two-dimensional electron 
gas, the electrons follow cyclotron orbits, and their allowed energies are quan- 
tized into values known as Landau levels [39.. At low magnetic field, the Lan- 
dau levels give rise to quantum oscillations including the de Haas- van Alphen 
effect and the Shubnikov-de Haas effect. At higher fields, the discrete Lan- 
dau level spectrum is manifest in the integer quantum Hall effect 40, 4T]l42]. 
a quantization of Hall conductivity into integer values of the quantum of 
conductivity e 2 /h. For monolayer graphene, the Landau level spectrum was 
calculated over fifty years ago by McClure [43] , and the integer quantum Hall 
effect was observed [15fT6] and studied theoretically [13,45,46,30,47] in recent 
years. The chiral nature of electrons in graphene results in an unusual sequenc- 
ing of the quantized plateaus of the Hall conductivity. In bilayer graphene, 
the experimental observation of the integer quantum Hall effect |17] and cal- 
culation of the Landau level spectrum [21) revealed further unusual features 
related to the chirality of electrons. 

1.8.1 The Landau level spectrum of monolayer graphene 

We consider a magnetic field perpendicular to the graphene sheet B = 
(0,0, —B) where B — |B|. The Dirac-like Hamiltonian of monolayer graphene 
(|1.52p may be written as 
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Hi,k+ = v 



vrt 

7T 



Hi. 



7T 
7ft 



{ 



Pa; " 
= Pr 



- 



(1.72) 



in the vicinity of corners of the Brillouin zone K+ and K _ , respectively. The 
off-diagonal elements of the Hamiltonian (|1.72[) contain operators 7r = p x + «Pj, 
and 7r' = p x — ipy, where, in the presence of a magnetic field, the operator 
p = (p x ,p y ) = —ihV + eA. Here A is the vector potential and the charge of 
the electron is — e. 

Using the Landau gauge A = (0, —Bx, 0) preserves translational invari- 
ance in the y direction, so that eigenstates may be written in terms of states 
that are plane waves in the y direction and harmonic oscillator states in the 
x direction [iTlH2"] . 



exp 



x 

As 



n 



: PyV 



(1.73) 



Here, T-Li are Hermite polynomials of order £, for integer £ > 0, and the 
normalization constant is Ag = 1/ \J2 1 £\^/tt. The magnetic length As, and a 
related energy scale r, arc defined as 



Ar = 



Ar 



(1.74) 



zBx and = 



With this choice of vector potential, it — —ifid x + hd y 
—ihd x — hdy + ieBx. Acting on the harmonic oscillator states (jl.73|) gives 



^4>i = 



i-i 



(1.75) 
(1.76) 



and ir<f>o — 0. These equations indicate that operators tt and are pro- 
portional to lowering and raising operators of the harmonic oscillator states 
<fie. The Landau level spectrum is, therefore, straightforward to calculate 
@3l|44,45 . At the first valley, K + , the Landau level energies and eigenstates 
of Hi k + arc 



K + ,£>1: 
K+,£ = 0: 



El± = ± 



E o =0, 



= -( ' 

V2 V=R<i 



, (1.77) 



(1.78) 



where ± refer to the conduction and valence bands, respectively. Equa- 
tion (|1.77p describes an electron (plus sign) and a hole (minus sign) series 
of energy levels, with prefactor r = ^/2Hv/\b (|1.74|) . proportional to the 
square root of the magnetic field. In addition, there is a special level (|1.78[) 
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fixed at zero energy that arises from the presence of the lowering operator in 
the Hamiltonian, 7r0 o = 0. The corresponding eigenfunction tp has non-zero 
amplitude on the A sublattice, but its amplitude is zero on the B sublattice. 
The form (|1.72p of the Hamiltonian Hx t K_ at the second valley, shows 
that its spectrum is degenerate with that at K + , with the role of the A and 
B sublattices reversed: 

K_,£>1: E u± 
K-,£ = 0: E a 

Thus, the eigenfunction -0o of the zero-energy level has zero amplitude on the 
B sublattice at valley K + and zero amplitude on the A sublattice at K- . If we 
take into account electronic spin, which contributes a twofold degeneracy of 
the energy levels, as well as valley degeneracy, then the Landau level spectrum 
of monolayer graphene consists of fourfold-degenerate Landau levels. 

1.8.2 The integer quantum Hall effect in monolayer graphene 

In this section, we describe how the Landau level spectrum of graphene is 
reflected in the dependence of the Hall conductivity a xy (n) on carrier density 
n. In conventional two-dimensional semiconductor systems, in the absence 
of any Berry's phase effects, the Landau level spectrum is given by Ei = 
hu) c (£ + 1/2), i > 0, where u c = eB/m is the cyclotron frequency [4Lll42| . 
Here, the lowest state lies at finite energy Eq = ?kj c /2. If the system has an 
additional degeneracy g (for example, g = 2 for spin), then plateaus |4QII41ll42j 
occur at quantized a xy values of N(ge 2 /h) where N is an integer and e 2 /h is 
the quantum value of conductance, i.e. each step between adjacent plateaus 
has height ge 2 /h, Fig. |1.12f a). Each a xy step coincides with the crossing of 
a Landau level on the density axis. Since the maximum carrier density per 
Landau level is gB/ipo, where ipa = h/e is the flux quantum, the distance 
between the a xy steps on the density axis is gB / ip^ . 

As described above, monolayer graphene has fourfold (spin and valley) 
degenerate Landau levels Et t ± — ±V2£?iv/Xb for £ > 1 and Eq = 0. The 
Hall conductivity (T xy (n), Fig. 11.12T b). displays a series of quantized plateaus 
separated by steps of size 4e 2 //i, as in the conventional case, but the plateaus 
occur at half-integer values of 4e 2 //i rather than integer ones: 

1 /4e 2 \ 
a xy = --(2N + l) ^— J , (1.81) 

where N is an integer, as observed experimentally [151116] and described theo- 
retically [13.45.46.30.47]. This unusual sequencing of a xy plateaus is explained 
by the presence of the fourfold-degenerate Landau level Eq fixed at zero en- 
ergy. Since it lies at the boundary between the electron and hole gases, it 



= ± 



= 0, 



1 ( ±i(f>£-i 
V2\ <t>e 



(1.79) 
(1.80) 
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(a) CT xy (geVh) (b) a xy (4e7h) (c) CT xy (4eVh) 




Fig. 1.12. Schematic representation of three types of integer quantum Hall effect, 
showing the density dependence of the Hall conductivity cr my (n): (a) conventional 
two-dimensional semiconductor systems with additional system degeneracy g; (b) 
monolayer graphene; (c) bilayer graphene. Here, B is the magnitude of the magnetic 
field and tpo = h/e is the flux quantum. 



creates a step in a xy of 4e 2 //i at zero density. Each Landau level in monolayer 
graphene is fourfold degenerate, including the zero energy one, so the distance 
between each a xy step on the density axis is AB/ipo, i.e. the steps occur at 
densities equal to integer values of AB/tpQ. 

1.8.3 The Landau level spectrum of bilayer graphene 

In the presence of a perpendicular magnetic field, the Hamiltonian (|1 .69|) 
describing massive chiral electrons in bilayer graphene may be written as 

H ^=-L(^ 2 )> ((*?)' o)' 

in the vicinity of corners of the Brillouin zone K + and respectively. 
Using the action of operators tt and 7P on the harmonic oscillator states 
<fie, (|1.75[) and (|1.76[) . the Landau level spectrum of bilayer graphene may 
be calculated [5T]. At the first valley, K + , the Landau level energies and 
eigenstates of H 2y K + are 

h 2 

K+,l=l: Ei=0, 

K+,£ = 0: £ = 0, 

where ± refer to the conduction and valence bands, respectively. Equa- 
tion (jl.83|) describes an electron (plus sign) and a hole (minus sign) series of 



^o=( o °), (1-85) 
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energy levels. The prefactor h 2 /(mX 2 B ) is proportional to the magnetic field, 
and it may equivalently be written as r 2 /ji or as hui c where w c = eB/m. For 
high levels, £ 3> 1, the spectrum consists of approximately equidistant lev- 
els with spacing Huj c . Note, however, that we are considering the low-energy 
Hamiltonian, so that the above spectrum is only valid for sufficiently small 
level index and magnetic field yflr <C 71. As well as the field-dependent levels, 
there are two special levels, (|1.84|) and (|1.85|) . fixed at zero energy. There are 
two zero-energy levels because of the presence of the square of the lowering 
operator in the Hamiltonian. It may act on the oscillator ground state to give 
zero energy, 7r 2 o = 0, (II -85[) . but also on the first excited state to give zero 
energy, 7r 2 0i = 0, (|1.84l) . The corresponding eigenfunctions ipo and %pi have 
non-zero amplitude on the Al sublattice, that lies on the bottom layer, but 
their amplitude is zero on the B2 sublattice. 

The form (| 1 . 82[) of the Hamiltonian i?2,if_ at the second valley, K_, shows 
that its spectrum is degenerate with that at K + with the role of the Al and 
B2 sublattices reversed. It may be expressed as i?2,if_ = &xH2,k+Cx so that 
ipt,±( K -) = a xil>t,±(K+), tpi(K-) = a x ipi(K + ), and ipo(K-) = a x ip (K+). 
Thus, the eigenfunctions tpo and ip\ of the zero-energy levels have zero am- 
plitude on the B2 sublattice at valley K + and zero amplitude on the Al 
sublattice at If we take into account electronic spin, which contributes 
a twofold degeneracy of the energy levels, as well as valley degeneracy, then 
the Landau level spectrum of bilayer graphene consists of fourfold degenerate 
Landau levels, except for the zero-energy levels which are eightfold degener- 
ate. This doubling of the degeneracy of the zero-energy levels is reflected in 
the density dependence of the Hall conductivity. 

1.8.4 The integer quantum Hall effect in bilayer graphene 

The Hall conductivity a xy (n) of bilayer graphene, Fig. ll.12T c). displays a series 
of quantized plateaus occurring at integer values of 4e 2 /h that is practically 
the same as in the conventional case, Fig. 11.12( a). with degeneracy per level 
g = 4 accounting for spin and valleys. However, there is a step of size 8e 2 /h in 
a xy across zero density in bilayer graphene |17U21j . This unusual behavior is 
explained by the eightfold degeneracy of the zero-energy Landau levels. Their 
presence creates a step in a xy at zero density, as in monolayer graphene, but 
owing to the doubled degeneracy as compared to other levels, it requires twice 
as many carriers to fill them. Thus, the transition between the corresponding 
plateaus is twice as wide in density, 8B/ipo as compared to AB/ipo, and the 
step in a xy between the plateaus must be twice as high, 8e 2 /h instead of 
4e 2 /h. This demonstrates that, although Berry's phase 27r is not reflected in 
the sequencing of quantum Hall plateaus at high density, it has a consequence 
in the quantum limit of zero density, as observed experimentally |17j . 

Here, we showed that the chiral Hamiltonians of monolayer and bilayer 
graphene corresponding to Berry's phase n and 2-7T, respectively, have asso- 
ciated four- and eight-fold degenerate zero-energy Landau levels, producing 
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steps of four and eight times the conductance quantum e 2 /h in the Hall con- 
ductivity across zero density [TMTfilfTT] . In our discussion, we neglected inter- 
action effects and we assumed that any valley and spin splitting, or splitting 
of the £ = 1 and £ = levels in bilayer graphene, are negligible as compared 
to temperature and level broadening. 



1.9 Trigonal warping in graphene 

So far, we have described the tight-binding model of graphene and showed 
that the low-energy Hamiltonians of monolayer and bilayer graphene sup- 
port chiral electrons with unusual properties. There are, however, additional 
contributions to the Hamiltonians that perturb this simple picture. In this 
section, we focus on one of them, known in the graphite literature as trigonal 
warping [H[53[M1[2S1I1] • 



1.9.1 Trigonal warping in monolayer graphene 

The band structure of monolayer graphene, shown in Fig. ll.4l is approximately 
linear in the vicinity of zero energy, but it shows deviations away from linear 
behavior at higher energy. In deriving the Dirac-like Hamiltonian of monolayer 
graphene (jl.52j) . we kept only linear terms in the momentum p = Kk — TiK.^ 
measured with respect to the point. If we retain quadratic terms in p, then 
the function /(k), (jl.51[) . describing coupling between the A and B sublattices 
becomes 



V3a , x a 



i 

2 



/ (k) « -^-(&>* - Wy) + ^2 + iPv) > (i- 86 ) 

where pa/h <C 1. Using this approximate expression, the Dirac-like Hamilto- 
nian (11.521) in the vicinity of point is modified [5] as 

Hl ' S - V Up X + iPy ) M Ufe-^) 2 J' (L87) 

where parameter /i = 7oa 2 /(8?i 2 ). The corresponding energy eigenvalues are 

E± = ±^v 2 p 2 ~ 2^vp 3 cos Zip + /xV , (1.88) 

For small momentum near the K point, pa/h <C 1, the terms containing 
parameter /i are a small perturbation because /jp 2 /(vp) — pa/(4y3h). They 
contribute to a weak triangular deformation of the Fermi circle that becomes 
stronger as the momentum p becomes larger. Figure fl . 131 shows the trigonal 
warping of the Fermi circle near point K+, obtained by plotting (|1.88|) for 
constant energy E = 0.57o- The presence of the valley index £ = ±1 in the 
angular term of (| 1 .88[) means that the orientation of the trigonal warping at 
the second valley K_ is reversed. 
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Fig. 1.13. Trigonal warping in monolayer graphene. The solid line shows the isoen- 
ergetic line E = O.570 in the vicinity of the valley K+ using Eq. (JT755J , the dashed 
line shows the circular isoenergetic line obtained by neglecting trigonal warping 



1.9.2 Trigonal warping and Lifshitz transition in bilayer graphene 

In deriving the low-energy Hamiltonian of bilayer graphene (|1.69j) . the lin- 
ear approximation of /(k) (|1.51j) in the vicinity of the K point was used. 
Taking into account quadratic terms in /(k) would produce higher-order in 
momentum contributions to (|1.69[) . that would tend to be relevant at large 
momentum p. There is, however, an additional interlayer coupling in bilayer 
graphene that contributes to trigonal warping and tends to be relevant at 
small momentum p, i.e. at low energy and very close to the K point. 

The additional coupling is a skew interlayer coupling between p z orbitals 
on atomic sites Al and B2, Fig. 11.141 denoted 73. For each Al site, there are 
three B2 sites nearby. A calculation of the matrix element between Al and 
B2 sites in the tight-binding model proceeds in a similar way as that between 
adjacent A and B sites in monolayer graphene, as described in Sect. 11.4.21 
Then, the effective Hamiltonian in a basis with components Al, Bl, A2, B2, 
for the four low-energy bands of bilayer graphene (11.58)) is [3T] : 



(i = 0. 



H = 



( 

-7o/* (k) 



7o/ (k) 







7i 








73/* (k)\ 







(1.89) 



V -73/ (k) 







7l 




7o/* (k) 



7o/ (k) 
/ 



where 



73 = - {<Pai (r - Rai) |H|0s2 (r - Rbz)) 



(1.90) 
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Fig. 1.14. Schematic representation of the crystal structure of A_B-stacked bilayer 
graphene illustrating skew interlayer coupling 73 (dashed lines) between p z orbitals 
on sites Al (white) and B2 (grey) 



The 73 term is relevant at low energy because it is a direct coupling be- 
tween the Al and B2 orbitals that form the two low-energy bands. Thus, 
using the linear- in- momentum approximation (|1.51[) . terms such as 73/ (k) w 
—V3 {^p x ~ ipy) appear in the two-component Hamiltonian written in basis 
cai, cb2- Equation (|1.69[) is modified as [21] 



#2,£ = «3 







&x - iPy 



1 

2m 







{&x -iPy) 





(1.91) 



where W3 = \/3a73 / (2h) and m = r yi/(2v 2 ). Taking into account trigonal 
warping, the low-energy Hamiltonian of bilayer graphene (jl.9ip resembles that 
of monolayer graphene (II .8T[) . The principle difference lies in the magnitude of 
the parameters. Since 73 = 0.315 eV [4] is an order of magnitude less than 70 = 
3.033 eV |llj . then w 3 <C v. Thus, the linear term dominates in monolayers 
and the quadratic term dominates in bilayers over a broad range of energy. 
The energy eigenvalues of H 2 ^, (|1.91[) , are 



E ± =±dvip>-£ 



.v 3 p° 



cos 3f + — 
2m 



(1.92) 



for energies \E±\ <C 71. Over a range of energy, the term independent 
of V3 dominates, and the V3 dependent terms produce trigonal warping of 
the isoenergetic line in the vicinity of each K point. The effect of trigo- 
nal warping increases as the energy is lowered, until, at very low energies 
El ~ jji{ v 3/ v ) 2 ~ lmeV, it leads to a Lifshitz transition [48]: the isoen- 
ergetic line breaks into four parts [H [23 HI [23 H ED 32 [S^- There is one 
'central' part, centered on the K point (p = 0), that is approximately circular 
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with area A c ~ ttE 2 / (Tlv^) 2 . In addition, there are three 'leg' parts that are 
elliptical with area Ag ~ \A C . Each ellipse has its major axis separated by 
angle 2-7r/3 from the major axes of the other leg parts, as measured from the 
K point, with the ellipse centered on |p| = J1V3/V 2 . 

Here, we have described the low-energy band structure of monolayer and 
bilayer graphene within a simple tight-binding model, including a Lifshitz 
transition in bilayer graphene at very low energy El ~ 1 meV. It is quite 
possible that electron-electron interactions have a dramatic effect on the band 
structure of bilayer graphene, producing qualitatively different features at low 
energy [5T1IS1IS3I11IS3ISE] • 



1.10 Tuneable band gap in bilayer graphene 

1.10.1 Asymmetry gap in the band structure of bilayer graphene 

In graphene monolayers and bilayers, a combination of space and time in- 
version symmetry [57 guarantees the existence of a gapless band structure 
exactly at the K point, i.e. the A and B sublattices (Al and B2 in bilay- 
ers) are identical, leading to degeneracy of the states they support at the K 
point. Breaking inversion symmetry by, say, fixing the two sublattice sites 
to be at different energies, would lead to a gap between the conduction and 
valence bands at the K point. In monolayer graphene, breaking the A/B sub- 
lattice symmetry in a controllable way is very difficult: it would require a 
periodic potential because A and B are adjacent sites on the same layer. In 
bilayer graphene, however, the A\ and B2 sublattices lie on different layers 
and, thus, breaking the symmetry and opening a band gap may be achieved 
by doping or gating. Band-gap opening in bilayer graphene has recently been 
studied both theoretically [ITllMlESEHllMllSIllMllMlEaEl] and in a range of 
different experiments [l^lHlIMllSSIHSEilMllIi^ 

If we introduce an asymmetry parameter A = €2 — ei describing the dif- 
ference between on-site energies in the two layers, £,42 = e_B2 = £2 = ^A, 
€ai = cbi = ei = then the transfer integral matrix of bilayer graphene 

fT55]) . in a basis with components Al, Bl, A2, B2, becomes 



H = 



( -\A -70/ (k) \ 

-7o/*(k) -\A 71 
71 k A -7o/ (k) 

V -70/* (k) \A J 



(1.93) 



The band structure may be determined by solving the secular equation 
det (H — EjS) — using overlap matrix 5, (11.571) . It is plotted in Fig. 11.151 for 
parameter values 70 = 3.033 eV, so = 0.129 and A = 71 = 0.39 eV. A band 
gap appears between the conduction and valence bands near the K points 
(left inset in Fig. ll.15jl . 
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Fig. 1.15. The low-energy band structure of bilayer graphene in the presence of 
interlayer asymmetry A. Parameter values are 70 = 3.033 eV, so = 0.129, A = 71 = 
0.39 eV. The plot shows the bands calculated along the k x axis intersecting points 
r, and K+ in the Brillouin zone, shown as the dotted line in the right inset. 
The left inset shows the band structure in the vicinity of the point K- . 



To develop an analytic description of the bands at low energy, we neglect 
non-orthogonality of the orbitals on adjacent sites, so that the overlap matrix 
S, (|1.57[) , becomes a unit matrix. Then, the bands at low energy are described 
by Hamiltonian, (|1.93l) . with eigenvalues jH] given by 



E 



(a) 



± 




(1.94) 



where a = 1 for the split bands and a = — 1 for the low-energy bands. 
Here, we used the linear approximation / (k) m —v (£p x — ip y ) /70, (|1.51[) . so 
that 7o|/(k) I w vp. Eigenvalues E^ 1 ^ describe the low-energy bands split 
by a gap. They have a distinctive 'Mexican hat' shape, shown in the left 
inset in Fig. 11.151 The separation between the bands exactly at the K point, 
(p = 0) — E^'ip = 0), is equal to \A\, but the true value of the band 
gap A g occurs at non-zero value of the momentum p g away from the K point, 




E { S 1} ( Pg ) 



1! 



(1.95) 
(1.96) 
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For moderate values of the asymmetry parameter, \A\ <C 71, then the band 
gap A g ss \A\, but for extremely large values, \A\ >• jx, the gap saturates 
A g ss 7!, where 71 is of the order of three to four hundred meV. The value 
of the asymmetry parameter A and bandgap A g may be tuned using an 
external gate potential, but the ability of an external gate to induce a potential 
asymmetry between the layers of the bilayer depends on screening by the 
electrons in bilayer graphene, as discussed in the following. 

1.10.2 Self-consistent model of screening in bilayer graphene 
Introduction 

The influence of screening on band-gap opening in bilayer graphene has been 
modeled using the tight-binding model and Hartree theory [28 , 60,62, 65J[SS] , 
and this simple analytic model is in good qualitative agreement with density 
functional theory 60,64_ and experiments [62 67 63 69,71. Recently, the tight- 
binding model and Hartree theory approach has been applied to graphene 
trilayers and multilayers [75j [76l [77l [78] . Here, we review the tight-binding 
model and Hartree theory approach which provides analytical formulae that 
serve to illustrate the pertinent physics. We will use the SI system of units 
throughout, and adopt the convention that the charge of the electron is — e 
where the quantum of charge e > 0. 

Using elementary electrostatics, it is possible to relate the asymmetry pa- 
rameter A — 62 — t\ to the distribution of electronic density over the bilayer 
system in the presence of external gates, but the density itself depends explic- 
itly on A because of the effect A has on the band structure, (|1.94|) . Therefore, 
the problem requires a self-consistent calculation of density and A, leading to 
a determination of the gate-dependence of the gap A g . 

The model assumes that bilayer graphene consists of two parallel conduct- 
ing plates located at x = — c /2 and +co/2, where Co is the interlayer spacing, 
as illustrated in Fig. 11.161 The two layers support electron densities n\, 122, 
respectively, corresponding to charge densities o\ = —eni, 02 = —eri2, and 
the permittivity of the bilayer interlayer space is e r (neglecting the screening 
effect of 7r-band electrons that we explicitly take into account here). We con- 
sider the combined effect of a back and top gate, with the back (top) gate at 
x = —Lb (x = +L t ), held at potential V& (V t ), separated from the bilayer by 
a dielectric medium with relative permittivity {st). In addition, we include 
the influence of additional background charge near the bilayer with density 
ribo on the back-gate side and nto on the top-gate side, yielding charge den- 
sities (Tf,o — erif,o and &m = e^-to where nbo and rito are positive for positive 
charge. 

Electrostatics 

Applying Gauss's Law firstly to a Gaussian surface enclosing cross-sectional 
area A of both layers of the bilayer and, secondly, to a Gaussian surface 
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Fig. 1.16. Schematic of bilayer graphene in the presence of back and top gates. 
Bilayer graphene consists of two parallel conducting plates with respective electron 
densities m, ri2 located at x = — co/2 and +co/2, respectively, where Co is the 
interlayer spacing, and e r is the permittivity of the bilayer interlayer space. The 
back (top) gate at x — —Lb (x = +Lt), held at potential Vb (Vt), is separated from 
the bilayer by a dielectric medium with relative permittivity Eb (et). Dashed lines 
indicate additional background charge near the bilayer with charge densities Obo and 
(Jto on the back-gate and top-gate side, respectively. 



enclosing one layer only yields 

- e a e b E b A + e Q e t E t A = -e (m + n 2 - n ba - n m ) A , (1-97) 
-e e r EA + e e t E t A = -e (n 2 - n t0 ) A . (1.98) 

The electric fields may be related to potential differences, 

E b a V b /L b , Et pb -Vt/Lt , (1.99) 
En(Vi- V 3 ) /cq = A/(ec ) . (1.100) 

and, when substituted into (|1.97p and (|1.98[) . they give 

n = ni+n 2 = 1 h n b0 + n t0 , (1.101) 

eL b eL t 

A = -Ei^ eVt + ^ {n . 2 - nm) . (i.io2) 

E r Lit SoS r 

The first equation, (| 1 . 101[) , relates the total density of 7r-band electrons n — 
ni+n 2 on the bilayer to the gate potentials, generalizing the case of monolayer 
graphene Q3] . The second equation, (|1.102p , gives the value of the asymmetry 
parameter. Using (|1.101[) . it may be written in a slightly different way: 

^ = ^ ext + ^7i (n2 ~ ni) , (1-103) 

le b c 1 e t c (n b0 - n t0 ) 
^ext = 7: — i-eVb - - — — eVt + Aji , (1.104) 

Ie t L b I e r Lt n±_ 



where parameters n± and A are defined as 
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7i , c e 2 7i _ c e 2 n_L n ln -s 

= — ^ — , A — 5 = . (1.105) 

Trh v 2 2nh V 2 £ £r ^ll^r 

The first term in (|1.103|) is A ex t, the value of A if screening were negligible, 
as determined by a difference between the gate potentials, (|1.104jl . Equa- 
tions (|1. 10111. 104p show that the effect of the background densities ribo and 
nto m ay be absorbed in a shift of the gate potentials V& and Vt , respectively. 

The second term in (|1.103l) indicates the influence of screening by elec- 
trons on the bilayer where n± is the characteristic density scale and A is a 
dimensionless parameter indicating the strength of interlayer screening. Using 
7i = 0.39eV and v = 1.0 X 10 6 ms _1 gives nj_ = 1.1 x 10 13 cm -2 . For interlayer 
spacing c = 3. 35 A and dielectric constant e r ~ 1, then A ~ 1, indicating 
that screening is an important effect. 



Layer densities 

Equation (|1 .103[) uses electrostatics to relate A to the electronic densities n\ 
and ri2 on the individual layers. The second ingredient of the self-consistent 
analysis are expressions for m and ri2 in terms of A, taking into account the 
electronic band structure of bilayer graphene. The densities are determined 
by an integral with respect to momentum over the circular Fermi surface 

™K2) = "L fp d P DVUl(2) (P)\ 2 + I^Bl(2)b)| 2 ] (1-106) 

irn J 

where a factor of four is included to take into account spin and valley degener- 
acy. Using the four-component Hamiltonian (|1.93[) , with linear approximation 
/ (k) w — v (£p x — ip y ) /70, it is possible to determine the wave function am- 
plitudes on the four atomic sites [28] to find 



(1.107) 



(E 2 - A 2 1 A) T 2v 2 p 2 EA - v 4 p 4 



(E 2 - A 2 /Ay + v 2 p 2 A 2 - dV 



where the minus (plus) sign is for the first (second) layer and E is the band 
energy. 

For simplicity, we consider the Fermi level to lie within the lower con- 
duction band, but above the Mexican hat region, \A\/2 < Ep <C 71 • We 
approximate the dispersion relation, (ll.94[) . as ^ ss yj A 2 /A + v 4 p 4: /j 2 
which neglects features related to the Mexican hat. Then, the contribution to 
the layer densities from the partially-filled conduction band [2SH15] is given 
by 



,cb - n - n ^ A 
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where the total density n = p 2 F / irfi 2 . In addition, although the filled valence 
band doesn't contribute to a change in the total density n, it contributes 
towards the finite layer polarization in the presence of finite A which, to 
leading order in A, is given by 



vb 



Then, the total layer density, n\(2) 



-.cb 
H(2) 



471 

\A\ 



l l(2) 



is given by 



(1.109) 



M(2) 



n n±A , 
- =P In 

1 471 




ITT • (LUG) 



Self-consistent screening 

The density-dependence of the asymmetry parameter A and band gap A g 
are determined [281 65 by substituting the expression for the layer density, 
(fTTTOj) . into (|LT03| : 



A (n) « Z\ c 




(1.111) 



with Z\ ext given by (|1.104|) . The logarithmic term describes the influence of 
screening: when this term is much smaller than unity, screening is negligible 
and A « A ext , whereas when the logarithmic term is much larger than unity, 
screening is strong, \A\ <C |Z\ cx t|- The magnitude of the logarithmic term is 
proportional to the screening parameter A. As discussed earlier, A ~ 1 in 
bilayer graphene, so it is necessary to take account of the density dependence 
of the logarithmic term in (ll.lllj) . 

To understand the density dependence of A, let us consider bilayer 
graphene in the presence of a single back gate, Vt = — n t Q = 0. This 
is a common situation for experiments with exfoliated graphene on a silicon 
substrate [T4 t [T5 l H6 l ll7j. Then, the relation between density and gate voltage, 
(jl.lOip . becomes the same as in monolayer graphene [H], n — eo£bVb / (eLt,) . 
The expression for A ext , (|1.104[) . reduces to Z\ cxt = A'fin/n±, and the expres- 
sion for A (n) 1 (jl.llll) . simplifies [28] as 



A{n) 



Ajin 



2 \n± 



(1.112) 



The value of the true band gap A g (n) may be obtained using (11.96[) . A g = 
I-^ITi/v ^ 2 + 7i- Asymmetry parameter A(n) and band gap A g (n) are plot- 



ted in Fig. 11.171 as a function of density n. For large density, 



nj_, 



the 
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Fig. 1.17. Density-dependence of the band gap A g in bilayer graphene, in the 
presence of a single back gate: (a) asymmetry parameter A and gap A g = 
\A\>n/y/A 3 +f* for screening parameter A = 1; (b) band gap A g for different 
values of the screening parameter. Plots were made using Eqs. (|1.96f) and (|1.112l) . 

logarithmic term in (| 1 . 1 1 2[) is negligible and the asymmetry parameter is ap- 
proximately linear in density, A(n) ~ A'jin/n^. At low density, |n| <C n±, 
the logarithmic term is large, indicating that screening is strong, and the 
asymmetry parameter approaches A(n) ~ 2"/i(n/n±)/ hi(n±/\n\). The com- 
parison of A(n) and A g (n), Fig. ll.17T a). shows that, at low density |n| <§; n±, 
A g (n) w |^\(?t-)| and, asymptotically, A g (n) w 2^i{\n\/n^)/\n(n±/\n\). This 
is independent of the screening parameter A [see Fig. 11.17( b)]. The curves 
for different values of the screening parameter A, Fig. 11.17T b). illustrate that, 
even when \A\ is very large, \A\ 3> 71, A g saturates at the value of 71. 

In deriving the above expression for A(n), a number of approximations 
were made including simplifying the band structure [by omitting features re- 
lated to the Mexican hat or to other possible terms in the Hamiltonian (|1.93[) ]. 
neglecting screening due to other orbitals, and neglecting the effects of dis- 
order and electron-electron exchange and correlation. Nevertheless, it seems 
to be in good qualitative agreement with density functional theory calcula- 
tions [BTTllrJI] and experiments (see, for example, [62 ( 67, 63, 69, 71 ). 

1.11 Summary 

In this Chapter, some of the electronic properties of monolayer and bilayer 
graphene were described using the tight-binding model. Effective Hamiltoni- 
ans for low-energy electrons were derived, corresponding to massless chiral 
fermions in monolayers and massive chiral fermions in bilayers. Chirality in 
graphene is manifest in many electronic properties, including anisotropic scat- 
tering and an unusual sequence of plateaus in the quantum Hall effect. There 
are a number of additional contributions to the low-energy Hamiltonians of 
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graphene that influence chiral electrons and we focused on one of them, trig- 
onal warping, here. 

Comparison with experiments suggest that the tight-binding model gener- 
ally works very well in graphene. The model contains parameters, correspond- 
ing to the energies of atomic orbitals or to matrix elements describing hopping 
between atomic sites, that cannot be determined by the model. They must 
be estimated by an alternative theoretical method, such as density-functional 
theory, or they can be treated as fitting parameters to be determined by 
comparison with experiments. The simple model described in this Chapter is 
versatile and it serves as the starting point for a wide range of models en- 
capsulating advanced physical phenomena, including interaction effects, and 
the tight-binding model may be used to describe the electronic structure of 
multilayer graphene, too. Here, we described a different example: the use of 
the tight-binding model with Hartree theory to develop a simple model of 
screening by electrons in bilayer graphene in order to calculate the density 
dependence of the band gap induced by an external electric field. 
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